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We calculate the Mellin moments of next-to-next-to-leading order coefficient functions of the Drell-Yan and Higgs 
production cross sections. The results can be expressed in term of finite harmonic sums which are maximally 
threefold up to weight four. Various algebraic relations among these finite sums reduce the complexity of the 
results suitable for fast numerical evaluations. It is shown that only five non-trivial functions occur besides Euler's 
^-function in the representation of these Wilson coefficients. 



1. Introduction 

With the present pp-collider Tevatron and 
the upcoming large hadron collider (LHC) [5] at 
CERN, the need for precise predictions from the- 
ory have become more and more important. The 
Drell-Yan (DY) cross section which is known upto 
the NNLO level [3] not only tests the reliability 
of perturbative QCD (pQCD) but also reduces 
the uncertainties coming from theory in order to 
make background studies more reliable for new 
particle searches and physics beyond the Stan- 
dard Model. Similarly, higher order corrections 
to the total cross section 0] for the Higgs bo- 
son production make the predictions more reliable 
for phenomenological studies. In the following we 
will study the structure of these corrections in the 
Mellin space using various algebraic identities re- 
lating the resulting finite harmonic sums. This 
representation allows a considerable reduction of 
the set of basic functions needed to represent the 
2-loop Wilson coefficients. 



* Based on the talk given at 7th DESY Workshop on El- 
ementary Particle Theory, Loops and Legs in Quantum 
Field Theory, April 25 -30, 2004, Zinnowitz (Usedom Is- 
land), Germany 



2. Coefficient Functions 

Due to mass factorization, hadronic cross sec- 
tions such as for the DY-process and Higgs- 
boson production can be expressed in terms of 
Mellin convolutions of the perturbatively com- 
putable coefficient functions A a b(x, Q 2 , /J 2 ) and 
non-perturbative parton distributions, / a (x,/x 2 ), 
of incoming hadrons. 

a(x,Q 2 ) = [ — [ —f a ( Xl:fl 2 )f b ( X2:fl 2 ) 

J x XX Jj-L. x 2 

x A ab (^-,%). (1) 

Though the parton densities are not calculable 
in the pQCD, their evolution with respect to a 
scale is computable, thanks to the renormaliza- 
tion group (RG) equations governing mass fac- 
torization. Hence higher order corrections to 
hadronic reactions enter through two sources viz, 
the coefficient functions and the RG-equations of 
parton distribution functions. In this paper we 
will concentrate only on the coefficient functions. 
Details of the calculation will be given in 

The perturbatively calculable coefficient func- 
tions arc usually computed in terms of the scaling 
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variable x = Q 2 /s where s is the center of mass 
energy of the incoming partonic system. Q 2 is 
the invariant mass of the final di-lepton pair (for 
DY) or mass squared of the Higgs boson. The 
results are expressible in terms of polynomials in 
x, logarithms and Nielsen integrals |Hj defined by 

S n , P {*)= \ ' 1M ■ / -log-^log^l-^) (2) 
(n - iy.pl J Q z 

In practice, using (1), one performs the integra- 
tion over x\ and x 2 after folding the perturba- 
tively computed coefficient functions with the ap- 
propriate parton distributions. This involves fur- 
ther evaluation of various Nielsen integrals given 
in (2) or even more general functions. Due to this 
the complexity of the numerical evaluation of the 
hadronic cross sections is rather large. In the fol- 
lowing, we will present an alternative treatment 
of the evaluation of the total cross section upto 
NNLO level by working in Mellin space. Such 
techniques have been used in the past to com- 
pute deep inelastic scattering cross sections and 
indeed they are found to be most suitable for var- 
ious resummation programs [?] . Recently, similar 
work has also been performed for the Wilson co- 
efficients for polarized and unpolarized deeply in- 
elastic scattering |8I9) . 

3. Mellin Moments and Finite Harmonic 
Sums 

Using the Mellin transform ^U] of a given func- 
tion F{x) 



M[F]{N)= [ dxx N - 1 F(x) 
Jo 



(3) 



the cross section (1) in the Mellin space (iV-space) 
becomes 

M[a](N,Q 2 ) = M[f a ](N,n 2 )M[f b ](N,v 2 ) 



xM[A ab ](N,^) . (4) 



As is clear from the above equation, the convo- 
lutions of functions in x-space reduce to simple 
products of Mellin moments in iV-space. The 
Mellin moments of these functions can be ana- 
lytically continued JT] to complex values of N so 



that one can use various analyticity properties of 
these functions in the complex iV-plane to evalu- 
ate them efficiently. 

For our analysis, the starting point is eqn. (1) 
with given parton densities / a (x,/i 2 ) and known 
coefficient functions A ab (x, Q 2 ) computed upto 
NNLO in pQCD. We then compute the Mellin 
moments of these functions in iV-space and per- 
form their analytic continuation to complex val- 
ues of N. At the end, we substitute them back 
into eqn. (4) and perform the inverse Mellin trans- 
formation to arrive at the results in ir-space us- 
ing a suitable contour in the complex TV-plane. 
Compared to the direct numerical convolution in 
x-space the numerical computation using the TV- 
space results is much faster. 

Since we are dealing with the DY/Higgs total 
production cross sections, they depend only on 
two variables, the scaling variable x and the vir- 
tually Q 2 , the invariant mass of the final state. 
The set of functions contributing to the two-loop 
coefficient functions in x space contains about 80 
elements. Their representation in Mellin space 
has been given in m terms of a specific class 
of sums |12I13I14I15| . the finite harmonic sums. 

The finite harmonic sum of m-indices is defined 



AN) 



N n-i 

[sj g n(fci)] ni [sign(k 2 )] n2 
^ n^ 1 A nl fc2 ' 



Til —1 



no — 1 



E 



[sign(k m )] n " 



(•») 



with I, ki ^ 0. For example, 
'log 3 (l-z) 



M 



(1-x) 



(N) = -St(N-l) 



+ ls 2 (N-l)S 2 (N-l) + ^S 2 2 tX I) 



+25i(iV - 1)S 3 (N - 1) + -5 4 (iV - 1) . (6) 

Here the right hand side is a polynomial out of 
only single harmonic sums. Similarly, a combina- 
tion of various Nielsen integrals may have a simple 
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structure in iV-space ^2] : 



M 



i-x) + \ (2£i 2 (-x) log(l + x) + \og(x) 



xlog 2 (l + a;; 



(iV) = ^("^ C2lOS(2) + ^ C3 



-(-l)*(5i,_ a (JV) + ±6(Si(J\T) - 5_!(iV)) 
+ iC3-k 2 log(2))) (7) 



The single harmonic sums (m = 1) obey the fol- 
lowing integral representation 



N + ;(-Y) = 

£1 



(±s fc )* - 1 



(8) 



(-1) 



fe-i r 1 _ (±x) N - 1 

da; log _ (a;)- 



(fc-1)! J Q v "' X : 1 

Using the above representation for the single sums 
and the following formula, 

N 



(±x) 
fc=i 



dzlog 1 1 (z) " 



(9) 

one derives integral representations for the mul- 
tiple sums Sk 1 ...k m (N). For our purpose, we have 
made extensive use of the results given in |12| . 

4. Algebraic Relations 

Finite harmonic sums are related to each other 
by various algebraic relations |16I17I12| . These 
relations result from the shuffle product of two 
harmonic sums. Special cases follow from partial 
or complete index permutations. 

The permutation relation for sums with two in- 
dices is given by 

(N) + S n , m {N) = S m (N)S n (N) 

(AO , (10) 

with m An — sign{m}sign{n}(|m| + |n|). 

For harmonic sums with three indices the fol- 
lowing identity is used to derive various relations 



E 

perm{l,m,n} 



Si, m , n (N) = Sz(JV)S m (JV)S n (JV) 



AmAn (ii) 

inv perm{l,m,n} 

where "inv perm" means the invariant permuta- 
tions and "perm" means permutations. 

SiAiW - -2S 2 ,i,i(N) + S 3>1 (N) 

+Si(N)S 2 ,i(N) + S 2 , 2 (N) 

Si,i, 2 (N) = S 2 ,i,i(N) + ~(Si(N)(Si, 2 (N) 

-S 2A (N)) + S h3 (N) - S 3 ,i(N)) 

51, _ 2 ,i(iV) = -2£_ 2i i il (JV)+S_ 3 ,i(iV) 

+S 1 (A r )S-2,l(A f )+£-2, 2 (A r ) 

Si,i,- 2 (N) = S_ 2 ,i,i(W) + S- 2 (N)S 2 (N) 
S- 2 , 2 (N) - S- 2 (N)S 1A (N) 
+Si(N)S h - 2 (N) + Si,- 3 (N) 
-Si(N)S- 3 (N) (12) 

5. Coefficient Functions in A^-Space 

We have computed the Mellin moments for the 
DY and the Higgs coefficient functions using the 
list given in ^2] an d the algebraic identities given 
in the previous sections. For a general investiga- 
tion of these relations, see |14I18| . The algebraic 
relations considerably simply our final result. At 
the intermediate stages of the computation we en- 
counter various complicated sums such as 5i,_i,2, 

S-1,-1,-2, S_i ) _2,-l) S— 2,-1,-1, S 2 ,-i t i, Sl, 2 ,-l, 

5 2 , \-i, 5-1.1,2, S-i,i, 2 . At the end, most of these 
sums disappear leaving only few sums like £-2,1,1; 
£2,1,1 ■ We have listed below the sums that appear 
at the end of the computation. 

We obtain eight single sums S±i(N) with i = 
1...4, 



S- 4 (N) 
S- 3 (N) 



-1) 



il 



M 



log 3 (z) 



1 + a; 



TCI 
20 



-D n Im 



log 2 (a;) 



(JV + l)- i Ca 
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S- 2 (N) = (-1) N+1 M 

S-i(N) = (-1) N M 
S 4 {N) 

S a (N) 
S 3 (N) 
Si(N) 



log(x) 



1 + x 



(N + 1) - -C2 



1 + £ 



(TV + 1) - log(2) 



-71/ 



log 3 (^) 
1-x 

log 2 (^) 
1 - a; 



(N + 1) + H Ca 2 



(/V + 1) + Ca 



M 



log(^) 
1 — x 



(N + 1) + C 2 



= — M 



(7V + 1) 



(13) 



These harmonic sums can be solely expressed 
in terms of the Euler ^-function and the /3- 
function |17| and their derivatives, which is re- 
lated to the former combining two ^-functions 
with shifted argument. These functions repre- 
sent at the same time the analytic continuation 
of these harmonic sums. The following five dou- 
ble sums S_ 3 ,i(JV), 5_ 2 ,i(JV), S_ 2>2 (2V), S 2 ,i(/V), 
S 3 ,i(N) occur : 



$-3,1 (N) 



(—1) N M 



Ci 3 (x) 



1 + x 



(N + l) 



+C2S_2(/V)-C3S_l(iV)--CI 



+2£u(~) +^C3log(2) 
-I< 2 log 2 (2) + llog 4 (2) 



5_ 2ll (iV) = 



-(-1)"M 



Cl2{x) 



1 + x 



(N + l) 



5-2, 2 (/V) 



+C 2 log(2) 



-(— 1) JV M 



-io g (x)(/:i 2 ( a; ) + C2)) ( v - h 



+C2S_ 2 (/V)+2C35'-i(/V) 



40 s V2 



^ C 3log(2)+< 2 log 2 (2)-^ 



S 2 ,i(JV) 



M 



Ci2(x) 
1-x 



(JV + l) + C 2 5i(iV) 



£j 2 (x) log(x) 
1-x 



(7V+1) 



+< 2 S 2 (/V) - -S-(iV) - -S 4 (JV) 
20 S2 



(14) 



Furthermore, two triple sums S_ 2) i ) i(iV), 
52,1,1 (N) contribute : 

~Su(x) 



1 + x 



(iV + 1) 

i(-V)- ^'M^+gd 



+ ic3log(2) + ic 2 log 2 (2) 



1 



2^(2) 



5 2 ,i,i(iV) = M 



5i 2 (x) 



(N + l) 



1-x 

+CaSi(N) (15) 

From our final expression, we observe that only 
very few functions do finally contribute to the 
Wilson coefficients. We list them below for com- 
pleteness: 

n = 0, 1,2, d 

1 — x 1 + x 

£z 2 (x)log"(x) £z 2 (x)log n (x) 



1-x 



1 + x 



n = 0, 1 



5 



Si 2 (x) 



Ci 3 (x) 
1 + x 



Su(x) 



(16) 



Since functions weighted by a factor In™ (a;) are 
related to the un-weighted functions in Mellin 
space by 



M[\n n (x)f(x)](N) 



gn 



M[f(x)](N) 



(17) 



and single harmonic sums can be expressed by 
the well-known V'-functions only five non-trivial 
functions are needed to express the two-loop Wil- 
son coefficients for the polarized and unpolarized 
DY-process and the scattering cross sections of 
hadronic Higgs-boson and pseudoscalar Higgs- 
boson production. 

The next step is to multiply the appropriate 
parton densities computed in complex A-space 
with our coefficient functions obtained to invert 
back to izi-space for phenomenological studies. 

The Mellin inversion ^Uj of a function F(N) is 
given by 



F(x) = 



1 

2TTi 



C+200 



dNx~"F(N) 



(18) 



Here the parameter c is the intersection of the 
contour on the real axis and is chosen in such a 
way that the integral J dxx c F(x) is conver- 
gent. In other words, c should lie on the right of 
rightmost singularity of the function F(N). The 
shape of the contour can be deformed at our con- 
venience provided all singularities are covered. 

6. Conclusion 

We have systematically analyzed the mathemat- 
ical structure behind the NNLO coefficient func- 
tions for DY and Higgs production using Mellin 
moment techniques. Use of various algebraic 
identities which relate the finite harmonic sums 
in Mellin TV-space reduces the complexity of the 
results from around 80 functions to only five ba- 
sic functions, the ^-function and a few deriva- 
tives thereof. This is very useful not only to un- 
derstand the nature of higher order corrections 
but also to perform fast numerical calculations at 



high precision for phenomenological applications 
and fits to data. The same structures are found 
in the case of polarized and unpolarized 2-loop 
fragmentation functions 1201 • Together with the 
results of [8| it is now shown that these structures 
are in common for all known massless 2-loop Wil- 
son coefficients. 
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